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ABSTRACT 

We present results on the X-ray properties of IC 342 using a 10 ks XMM-Newton observation. Thirty-five 
sources are detected coincident with the disk of IC 342 (more than tripling the number known), of which ks 3 1 
are likely to be intrinsic to IC 342. This point-source population exhibits a diverse range of spectral properties 
and has an X-ray luminosity function slope and infrared luminosity comparable to that of starburst galaxies 
such as M82 and the Antennae, while its relative lack of extended X-ray emission is similar to the properties of 
quiescent spirals. Although we find no evidence for short-term variability in any of the X-ray sources, we do 
detect long-term variability between this observation and the 1991 ROSAT and 1993/2000 ASCA observations 
for five sources. Notably, the second most luminous source in IC 342 (X-2) is found to have an absorption- 
corrected 0.5-10 keV luminosity of 5.75 x 10 39 erg s" 1 and a spectrum best-fit by an absorbed multi-color 
accretion-disk model with T in = 2.17^ keV and Nn = (15.3+}|) x 10 21 cm . This is the lowest luminosity 
state observed for X-2 to date, although the slope of the spectrum is intermediate between the previously 
observed low/hard and high/soft states. IC 342 X-l, on the other hand, is found to be in an identical state to 
that observed in 2000 with ASCA. Assuming X-l is in an anomalous very high (VH) state, then either (1) X-l 
has remained in this state between 2000 and 2002, and is therefore the longest duration VH-state binary ever 
observed, or (2) it was simply caught in a VH state by chance in both the 2000 ASCA and 2002 XMM-Newton 
observations. We have also confirmed the ROSAT HRI result that the nucleus of IC 342 is made up of both 
point-like and extended emission, with the extended emission contributing k55% and «35% in the 0.3-2.0keV 
and 2.0-10.0 keV bands, respectively. The spectrum of the nucleus is best fit by an absorbed two-component 
model consisting of a thermal plasma with a temperature of kT ks 0.3 and a power law with a photon index of 
r k 2.53. The relative fluxes of the two spectral components suggest that the nucleus is complex, with a soft 
extended component contributing approximately half of the total luminosity. 

Subject headings: galaxies: active — galaxies: general — X-rays: binaries — X-rays: galaxies 



1. INTRODUCTION 

IC 342 is a nearby, nearly face-on (z ~ 20°) Scd galaxy with 
spiral arms undergoing moderate star formation and an in- 
tense starburst core. It has an absolute visual magnitude of 
My ~ -21.3, an optical D25 diameter of ks 29.7', and neutral 
hydrogen and total dynamical masses o f Mhi = 2.1 x 10 9 Mm 
andAf tnt = l.l x 10 11 M n , respectively (Tullv & Fishe rll988t 
ICrosthwaite et all 120001) . Although the Galactic absorption 
column de nsity toward IC 34 2 is relatively large (Nh = 3 x 
10 21 cm" 2 ; IStark etalJll992l) . it has been observed by Ein- 
stein, ROSAT, and ASCA. ROSAT High Resolution Imager 
(HRI) observations ($fwhm ~ 5") resolved the disk of the 
galaxy into a collection of ten point sources above a 0.1- 
2.5 keV flux limit of ks 3 x 10~ 14 erg cm" 2 s" 1 (or an 
absorption-corr ected luminosity limit of ks 2 x 10 37 erg s" 1 ; 
iBregman et aHH993l hereafter BCT93). 2 ROSAT and ASCA 
measurements indicate that the bulk of the X-ray emis- 
sion arises from four ultraluminous X-ray sources (ULXs, 
Lx ks 10 39 -10 40 erg s" 1 ) , one of which is coincident with 
the nucleus (BCT 93; lOkada etalJ 119981 hereafter 098; 
iKubota et al1l20Qll hereafter K01). Two ASCA observations 
(38 ks in 1993, 276 ks in 2000) have shown the ULXs to vary 
both i n intensity and spectral shape (098; K01; Sugiho et al. 
120011) . suggesting that they are single objects with masses of 
ks 10-100 M© (if radiating isotropically below their Edding- 

1 Department of Astronomy & Astrophysics, 525 Davey Lab, The Penn- 
sylvania State University, University Park, PA 16802. 

2 We adopt a distance of 3.3 ± 0.3 Mpc to IC 342, based on the period- 
luminosity relation of 20 Cepheid variables iSaha et al. 2002). Results from 
other papers have been corrected to this distance. 



ton limits). The nucleus, however, has remained constant and 
appears to be slightly extended (BCT93). Since this extended 
emission encompasses the surface-brightness enhancements 
seen at optical-to-radio wavelengths (ks 10-15") that are at- 
tributed to a high rate of ongoing star formation, BCT93 spec- 
ulated that the emission might be from a hot gas bubble. 

Here we report results from a 10 ks XMM-Newton obser- 
vation of IC 342. This observation provides the highest reso- 
lution and most sensitive hard X-ray imaging of IC 342 pub- 
lished to date, as well as a field of view which covers the entire 
extent of the galaxy. We describe the observation and reduc- 
tion of the X-ray data in Sj2] provide the basic X-ray properties 
of the detected sources in the context of previous results in £j3] 
and summarize our findings in 21 

2. XMM-NEWTON OBSERVATION AND DATA ANALYSIS 

IC 342 was observed on 2001 February 1 1 for «10 ks with 
the PN-CCD c amera JStruder et al l200ll and the two MOS- 
CCD cameras ( Tur ner et al J 1200 1 ) using the medium filter. 
The nucleus of IC 342 was placed at the aimpoint, allowing 
the entire optical extent of the galaxy to be imaged. The pro- 
cessing, screening, and analysis of the data were performed 
using the standard tools from SAS (v.5.3.3) and CIAO (v2.3), 
as well as custom IDL software. The raw MOS and PN data 
were initially processed using the standard epchain and em- 
chain pipeline scripts. Time intervals contaminated by soft- 
proton flares were identified using the background light curve 
in the 10-13 keV band. Approximately 200 s of exposure 
was excluded at the end of the observation due to flaring; the 
background for the rest of the observation remained relatively 
constant at < 0.2 counts s" 1 . Our final exposures were 9488 s, 
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FIG. 1 .— An XMM-Newton "false-color" 25' X 25' image of IC 342 with 
red, green, and blue representing 0.3-2.0 keV, 2.0-5.0 keV, and 5.0-10 keV 
emission, respectively. Prior to combination, each image was smoothed with 
an adaptive-kernel algorithm 1 Ebeling et al. 2003). Little diffuse emission 
is observed. In addition to a bright nuclear component (source 21), 34 off- 
nuclear point sources are detected. The sources have been labeled according 
to their source number in TablelTI 
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FIG. 2.— A 25' X 25' Digitized Sky Survey By-band image with XMM- 
Newton contours superimposed. Contours are calculated from the com- 
bined MOS image for the 0.3-10 keV band. The lowest contour indicates 
0.33 counts pixel - ' (Pb 3x the average background value), with each subse- 
quent contour increasing by a factor of four. All sources lie within the white 
12' radius circle shown. 



9507 s, and 4888 s for the MOS1, MOS2, and PN detectors, 
respectively. We selected only good event patterns: <12 for 
MOS imaging and spectroscopy, <12 for PN imaging, and 
<4 for PN spectroscopy. 

Source detection was initially performed on the MOS and 
PN images in the 0.3-10 keV band using the standard SAS 
eboxdetect and emldetect algorithms. A likelihood value of 10 
was imposed, corresponding to a significance level of «3.6<7. 
Thirty sources were detected, although we found that sev- 
eral apparent sources — including two with >100 MOS+PN 
counts — were detected by eboxdetect but rejected by emlde- 
tect? As a cross-check, we compared the sources found using 
eboxdete ct with those found using the CI AO wavdetect al- 
gorithm (Freeman et al. 2002), run on the merged MOS+PN 
image with a significance threshold of 1 x 10~ 7 . We found 
good agreement between the algorithms, and we have cho- 
sen the 35 sources common to both methods as our X-ray 
point-source sample. Table ^ lists the basic X-ray prop- 
erties of these sources, while Figure [2 indicates their po- 
sition in the X-ray image. Background-subtracted counts 
were extracted using circular apertures ranging in radius from 
w 17-68", depending on off-axis angle and the relative bright- 
ness of the source compared to the local background. Local 
backgrounds were determined after removal of the 35 point 
sources. Our background-subtracted detection limit is w 25 
MOS+PN counts, corresponding to an absorbed 0.3-10 keV 
flux of w (1-2) x 10" 14 erg cm" 2 s _1 . 

We have more than tripled the number of known X-ray 
sources spatially coincident with IC 342, compared to the 
ROSAT HRI observation reported by BCT93 of similar soft- 
band depth, thanks to the high-energy response of XMM- 
Newton. Figure |2 shows X-ray contours overlaid on an opti- 

3 These sources are noted in TablefTI 



cal image of IC 342, indicating that the majority of the X-ray 
sources coincide with the spiral arms of the galaxy. 

The nom inal astrometric a ccuracy for XMM-Newton is 
~ 4" (e.g., Jansen et al. 2001). We have improved upon this 
somewhat by matchi ng X-ray sources to optical sources from 
the r-band image of Sah a et aD d2002) and the United States 
Naval Observatory B1.0 catalog (USNOB1; iMonet et all 
120031) . 4 X-ray matching to optical sources was performed 
using a 4" matching radius, with 13 sources having poten- 
tial optical counterparts: the nucleus of IC 342, four spa- 
tially extended sources (potential Hll regions), five faint point 
sources with R ~ 21-25 within the spiral arms of IC 342, and 
three bright foreground stars with R ~ 11-13. Excluding the 
five faint point sources that potentially could be chance align- 
ments, we find average XMM-Newton Right Ascension and 
Declination offsets of w —3 ."4 and 0."4, respectively, result- 
ing in a registration accuracy of sal. "3 rms; positions have 
been corrected for these offsets throughout. These offsets do 
not change significantly (only by 0"3) if we align the X-ray 
and optical images using only the three bright stars and the 
nucleus (i.e., the most secure matches). Although we have 
limited statistics, we find no obvious rotation of the X-ray im- 
age to within our registration errors. 

Because IC 342 has a large angular extent and lies at low 
Galactic latitude (b = 10.°58), the chance superposition of 
Galactic X-ray sources and background AGN is a concern. 
Fortunately, IC 342 has a Galactic longitude of I = 138.°17, 
so our line-of-sight is relatively far from the Galactic Cen- 
ter and therefore intersects only a relatively small fraction of 
the Galaxy. Based on the numbers of serendipitous sources 
detected in four archival observations within 15° of IC 342 

4 The r-band image was astrometrically registered to an accuracy of f»0" 1 
RMS using 216 stars from the Tycho-2 and Guide Star v2.2 catalogs. 
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Table 1, X-ray Sources Spatially Coincident with IC 342 
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NOTE. — Column 1: ID number. Sources denoted by a "*" are not thought to be intrinsic to IC 342. Column 2: Source name given as XMMU JHHMMSS.S+DDMMSS. 
Sources denoted by a f were rejected by emldetect for no obvious reason (see S0. Column 3: Other names. R# and A# indicate the ROSAT and ASCA source numbers assigned by 
BCT93 and 098, respectively. Column 4: Off-axis angle in units of arcminutes. Column 5: Background-subtracted 0.3-10 keV counts, summed over all of the EPIC instruments 
(pn+MOSl+MOS2). Aperture photometry was performed using extraction radii ranging in radius from ~17-68", depending on off-axis angle and the relative brightness of the source 
compared to the local background (see 12*1 . The errors for th e source and background counts were computed following the method of Gehrels i 1 986 ) and were then combined following 
the "numerical method" described in §1.7.3 o flLvonsiH99ill . Column 6: Band ratios, defined as the ratio of counts between the hard and soft bands. The quoted band ratios have been 
corrected for differential vignetting between the hard band and soft band using the appropriate exposure maps. Errors for this quantity are calculated following the "numerical method" 
described in §1.7.3 of Lvons > 1991 ). Column 7: Observed, aperture -corrected 0.5-10 keV fluxes in units of 10~ 14 erg cm -2 s~ ] from the best-fit models to the XMM-Newton spectra (see 
Tabled- For faint sources (< 150 counts), fluxes were calculated assuming an average power-law spectrum with N\\ = 3.4 X 10 21 cm 2 and a photon index of F = 1.63 as determined 
from the stacked spectrum of the faint sourc es. Column 8: Ab sorption-corrected 0.5-10 keV luminosities in units of 10 38 erg s -1 from the best-fit models to the XMM-Newton spectra. 
Column 9: R magnitude as derived from the Saha et al. i 2002) r image (indicated by an "S") or USNOB1. The magnitude lower limits were estimated based on the faintest detectable 
sources within the vicinity of the X-ray source. Column 10: X-ray-to-optical flux ratio using columns 7 and 9 such that log(Fx/Fo)=log(/ r x)+5.5+W/2.5. Column 11: Observed, 
aperture -corrected 0.5-2 keV fluxes in units of 10~ 14 erg cm -2 s -1 from the ROSAT HRI derived using XIMAGE and PIMMS. The adopted spectral models were identical to those 
assumed in Column 5. Column 13: Possible optical ID and additional comments. Saha and USNOB1 offsets are provided in parentheses for sources with optical counterparts. 



(Proposal IDs 00011, 10946, 11220, 11229), and from ex- 
trapolation of the m edium-deep XMM-Newton log/V-logS of 
BaldFetal] J20021 i.e., « 60-100 sources per deg 2 down 
to a limiting 0.5-10 keV flux of « 1 x 10" 14 erg cm" 2 s" 1 ), 
we would expect ~ 7-12 foreground or background sources 
within a 12' radius of IC 342. Given that we are able to iden- 
tify four of our 35 sources with foreground stars (sources 2, 
11, 15, and 20), we suspect that the majority of the remain- 
ing 3 1 X-ray sources are associated with IC 342 based on their 
large X-ray-to-optical flux ratios; half have log(Fx/^b) ^ 0.5, 
and all but the nucleus have log(Fx/Fo) > -0.5. A com- 
parison of the X-ray-to-optical flux rat ios of these sources 
with Figure 1 of Maccac aro et alJ £1988, after accounting for 



average X-ray and optical band differences) suggests these 
sources are inconsistent with X-ray emission from all nor- 
mal stars and even some AGN (the most likely foreground 
and background contaminants, respectively), but they are 
fully consistent with X-ray binaries in IC 342. Additionally, 
sources that are significantly brighter than our limiting X-ray 
flux, and are therefore rarer, are even more likely to be associ- 
ated with IC 342. For instance, there is a < 0.1% chance that 
any of the four sources brighter than ss 5 x 10~ 13 erg cm -2 s -1 
are foreground or background X-ray sources. 

Point-source spectra were extracted using the apertures de- 
scribed above. Event PI values and photon energies were 
determined using the latest gain files appropriate for the ob- 
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servatio n. Th e X-ray spectra were analyzed using XSPEC 
iArnaudH 199 61 Unless stated otherwise, spectral parameter 
errors are for 90% confidence assuming one parameter of in- 
terest (Ax 2 = 2.7). The X-ray fluxes and absorption-corrected 
luminosities for all 35 sources were calculated from spectral 
fitting using XSPEC. None of the sources was affected by pile- 
up. For faint sources (< 150 counts) thought to be intrinsic to 
IC 342 (i.e., not obvious foreground stars), fluxes were calcu- 
lated assuming an average absorbed power-law spectrum with 
Mi = (3.4^)x 10 21 cm" 2 and a photon index of T= 1.63^;}|, 
as determined from the stacked spectrum of the faint sources. 
Fluxes for the two faint foreground stars (sources 2 and 15) 
were calculated assuming the spectrum of source 1 1 (see Ta- 
ble|2}- Note that if a faint source's spectrum deviates substan- 
tially from these average values, then the flux and absorption- 
corrected luminosity of the source may change as well. 

3. IC 342 X-RAY PROPERTIES 

The 3 1 point sources plausibly associated with IC 342 span 
over two orders of magnitude in X-ray luminosity and dis- 
play a diverse range of spectral properties, as implied from the 
band ratios listed in Table[T](B/? - 0. 1-3.4). The 0.5-10 keV 
X-ray luminosity function (XLF) of the sources is relatively 
flat down to Lx ~ 4 x 10 37 erg s" 1 with a slope of 0.5 ±0.1, as 
determined from maximum likelihood fitting of the differen- 
tial XLF. At luminosities fainter than this limit, incomplete- 
ness due to the sensitivity limit and absorption bias against 
soft sources makes the XLF slope uncertain. Somewhat sur- 
prisingly, we detect little diffuse X-ray e miss ion (only a com- 
pact region around the nucleus; see £|3.1> . In estimating 
the amount of diffuse emission, we only used the combined 
MOS1+MOS2 0.3-2.0 keV image for simplicity. We ex- 
tracted a total of 124 ±48 counts within a 5' radius aper- 
ture around the nucleus, using a vignetting-corrected 5-10' 
background annulus and masking out point sources using 
their 99% encircled-energy radius apertures. Assuming this 
emission would have an absorbed thermal plasma spectrum 
with a temperature of 0.5 ke V (1 keV) and column density 
of Nh = 6 x 10 21 cm" 2 (see fl3.3> . our 3er upper limit corre- 
sponds to a flux density of < 8.4 x 10~ 15 erg cm" 2 s" 1 arcmin" 2 
(< 9.9 x 10~ 16 erg cm" 2 s" 1 arcmin" 2 ) and an absorption- 
corrected 0.5-2.0 keV luminosity of Lx < 5.6 x 10 39 erg s" 1 
(L x < 3.0 x 10 38 ergs" 1 ). 

In addition to these global properties, the XMM-Newton 
observation provides constraints on the spatial, spectral, and 
temporal properties of the X-ray point-source population as 
detailed below. 

3.1. Spatial Extent of the Nuclear Source 
XMMU J034648.8+680546 (IC 342 X-3) 

This XMM-Newton observation offers an order of magni- 
tude improvement compared to the high-resolution ROSAT 
HRJ observations in the number of counts collected from the 
nucleus of IC 342, allowing us to improve measurements of 
the nuclear spatial extent reported by BCT93. The nuclear 
source appears to be azimuthally symmetric in the raw im- 
ages out to ~40", so we extracted soft and hard-band counts 
from it in l."l (1 MOS pixel) annular bins using the com- 
bined MOS1+MOS2 image; the spatial resolution of the PN 
is > 1.5x worse (XMM-Newton Users' Handbook §3.2.1.1) 
and is therefore not very useful for this purpose. We also gen- 
erated images of the MOS1+MOS2 PSF using the SAS cali- 
bration tool calview. Counts from these images were then 
extracted in an identical manner to the data and normalized by 
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FIG. 3.— The soft and hard-band MOS1+MOS2 radial profiles (small 
filled circles) of the IC 342 nucleus, compared to our simple model (solid 
curves) and the XMM-Newto n PS F (dashed curves) calculated at 1.5 keV 
and 6 keV, respectively; see i|3.1l for details. The average background has 
been added to both the model and the PSF. The shaded regions indicate the 
lcr deviation of the measured profile. Deviations from the PSF, particularly 
in the soft band, can be seen out to radii of ss20 and f»10 pixels, at which 
point the 0.07 counts pixel soft-band and 0.04 counts pixel -1 hard-band 
backgrounds begin to dominate (dotted lines), respectively. A MOS pixel is 
equivalent to I'.' I. 



the innermost nuclear annulus. The resulting radial surface- 
brightness distributions and la errors for the data are shown 
in Figure|3] 

The surface-brightness distribution of the nucleus lies 
above the background out to radii of ~ 15-20 pixels 
(~ 16-22") and is clearly extended compared to the PSF 
model. The accuracy of the MOS PSF model from calview 
was verified using source 26 (1' from the nucleus); devia- 
tions between this source and the PSF were much smaller than 
those seen in Figure[3]and were consistent with source statis- 
tics. The fact that the fraction of extended emission is larger 
in the soft band, even though the imaging quality is similar 
in both bands, further supports the reality of the extent. An 
emission model consisting of both a point source and a uni- 
form disk with an « 8" (128 pc) radius convolved with the 
instrument PSF represents the surface-brightness distribution 
well, although there is clearly still some residual scatter due 
to possible clumping or patchy absorption. This is especially 
noticeable in the hard band where individual point sources 
are more likely to dominate the emission. The relative con- 
tributions from the extended disk component are w 55% and 
w 35% in the soft and hard bands, respectively, in crude agree- 
ment with the BCT93 soft-band values (9" radius, 44%). The 
relative softness of the extended emission is consistent with 
its ex pected physical origin as supernova-heated hot gas (see 
also ©33 for X -ray spectral details). Moreover, the extent 
of the X-ray disk emission mirrors the surface brightness en- 
hancements seen at optical-to-radio wavelen gths (» 10-15"; 



opt 

e.g.JBoker et alJl997tlButa & McCallll999HSchinnerer et alJ 
2003) that are attributed to a starburst ring, further strength- 
ening our conclusions. 
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3.2. X-ray Timing Analysis 

3.2.1. Short-Term Variability 

Although the XMM-Newton observation of IC 342 was 
short, we do have adequate statistics to evaluate the short- 
term timing characteristics of some of the off-nuclear sources. 
To examine objectively the existence of any significant vari- 
ations in the count rate, we used the Kolmogorov-Smirnov 
(KS) statistic on the unbinned data to test the null hypothesis 
that each source plus background rate was constant over the 
duration of the exposure. None of the sources varied at the 
> 90% confidence level. 

3.2.2. Long-Term Variability 

Comparison of the XMM-Newton observation with the 
19 ks ROSAT HRI observation on 1991 February 13-16 al- 
lows us to constrain the soft-band long-term variability of the 
entire IC 342 X-ray source population, while comparison with 
the 38 ks ASCA observation on 1993 September 19 and the 
276 ks ASCA observation on 2000 February 24-March 1 al- 
lows us to constrain the full-band variability for X-l , X-2 , and 
X-3. The hard-band variability will be discussed in fl3.3l in the 
context of our X-ray spectral analyses. 

To measure the soft-band variability, we extracted the HRI 
counts at the position of each source using the XI MAGE 
software package, which corrects for vignetting and aper- 
ture losses. For detections, we used the values reported in 
BCT93, and for 3a upper limits we used the sosta algo- 
rithm. The HRI count rates were converted to 0.5-2 keV 
fluxes using PIMMS (column 10 of Table[l). The PIMMS in- 
put spectral model was set to the best- fit model found from 
our XSPEC spectral analyses in £13 . 3 1 We found that five 
sources were variable in the soft band over the 10-year period: 
XMMU J034556.0+680455 and XMMU J034646.0+680946 
faded by a factor of ~6 and ~3.3. compared to the ROSAT 
observation, respectively; XMMU J034719.2+681 129 and 
XMMU J034723.3+680857 brightened by at least factors 
of 1.5 and 1.4 compared to their ROSAT upper limits, re- 
spectively; and finally J03454 1.4+680241 (ROSAT source 2; 
BCT93) was not detected by our XMM-Newton observation at 
all, indicating a decrease of at least a factor of 5. Given the 
poor HRI statistics for this latter source, however, it may have 
been spurious, although strong variability or extreme X-ray 
spectral properties (i.e., a very soft spectrum) cannot be ruled 
out. 

3.3. X-ray Spectral Analysis 

For each source, we extracted the spectra and response ma- 
trices for each CCD camera separately and fit all three spectra 
with the same model in XSPEC. We initially fitted the spectra 
with absorbed power-law models. 5 The spectral properties of 
the 23 sources with fewer than 150 full -band counts are not 
well constrained and have been fitted only with a mean spec- 
trum as given in 5j2]to derive X-ray fluxes and luminosities. 
The 12 sources with more than 150 full -band counts all had 
adequate photon statistics for individual fits, although > 1000 
counts were typically necessary to rule out competing spectral 
models clearly. The best-fit spectral parameters for these 12 
sources are given in Table |2] All of the sources thought to be 
intrinsic to IC 342 exhibited spectral cutoffs below ~ 1 keV 
that are best fit with column densities equal to or larger than 

5 In most cases, a mekal thermal plasma model iMewe et al. 1985) was 
equally acceptable. 



the Galactic value (see fJD, consistent with these off-nuclear 
X-ray sources being located in star-forming regions and spi- 
ral arms within IC 342. We discuss our results for the four 
brightest X-ray sources in IC 342 with > 1000 counts below. 

3.3.1. XMMU J0345 56.0+68045 5 (IC 342 X-l ) 

IC 342 X-l is the brightest source in IC 342 and was 
recentl y found to b e coincident with a supernova remnant 
jRoberts et allEool . X-l has been extensively studied by 
ASCA, as it is considered to be one of the best examples 
of a massive (> 10M Q ) black-hole binary. 098 found this 
source in a luminous high state in 1993, with a spectrum best- 
fitted by an absorbe d multi-co lor disk model (diskbb model 
in XSPEC: lMitsuda et alJl 984) arising from an optically thick 
standard accretion disk around a black hole. The best-fit ab- 
sorbed diskbb parameters were Nh = (4.7 ± 0.3) x 10 21 cm" 2 
and r in = 1.77 ±0.05 keV, resulting in a 0.5-10 keV flux 
of Fo.5-iOkeV w 1.0 X 10~ n erg cm" 2 s . When it was 
revisited in 2000, K01 found the flux (Fo.s-iokev ~ 3.1 x 
10~ 12 erg cm" 2 s" 1 ) and spectrum of this source had changed 
dramatically, and it was best fit by an absorbed power-law 
model with Na = (6.4 ±0.7) x 10 21 cm" 2 and T = 1.73 ±0.06. 
An ionized Fe-K edge at 8.4 ± 0.3 keV with an optical depth 
of 0.9 ±0.5 was also seen. K01 attributed this spectral change 
to be a transition between a high/soft state and a low/hard 
state, as is observed i n many Galactic an d Magellanic black- 
hole binary systems. Kubota et al. (2002, hereafter K02) re- 
examined the 2000 ASCA data, finding that the spectrum de- 
viates from a power-law shape such that there appears to be 
significant softening above 5 keV. They proposed that the 
spectrum was not indicative of a low-hard state, but rather 
of an anomalous very high (VH) state as observed in some 
Galactic black-hole binaries, which was adequately fit by 
a strongly Comptonized optically thick accretion disk with 
7i„ = 1 .1 ±0.3 keV and = 2.2 ± 0.4. We note that the ASCA 
w2.'5 radius beam used to study X-l (compare with Figure 1 
of 098) would have also included sources 6, 9, and 1 1 from 
Tabled contaminating the ASCA source by ss 19% based on 
the relative XMM-Newton source fluxes. 

To compare with the ASCA spectral fits, we fitted both an 
absorbed diskbb model and an absorbed power law to the X- 
ray spectrum. The best-fit parameters for the absorbed diskbb 
model were T in = 1 .93^ keV and N u = (3.5^j) x 10 21 cm" 2 
(\ 2 = 206.9 for 193 degrees of freedom), while those for the 
absorbed power law were T = 1.68^'yg and Nh = (6.0^ 5) x 
10 21 cm" 2 (x 2 = 166.0 for 193 degrees of freedom). Al- 
though the diskbb fit is formally acceptable, there are system- 
atic residuals at the softest and hardest energies, and hence 
the power-law model appears to offer a superior fit to the data 
(see the X-l residuals in Figure 0}. The flux of the source is 
found to be Fo.5-10 keV ~ 2.15 x 10" 12 erg cm" 2 s" 1 The best- 
fit parameters from the power-law model are consistent with 
those obtained during the 2000 ASCA observation, suggesting 
that little changed spectrally in the «27.5 months between the 
latest ASCA observation and our XMM-Newton observation. 
The XMM-Newton flux, on the other hand, is as 30% lower 
than the 2000 ASCA measurement. This is a bit larger than the 
source confusion estimate made above (accounting for slight 
bandpass differences), but the difference is only marginally 
significant. Thus we conclude that X-l has remained rela- 
tively constant over the « 2 yrs between the latest ASCA and 
XMM-Newton observations and is therefore likely to be in a 
VH state similar to the one reported by K02. 
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Table 2. X-ray Spectral Analysis Results 



(1) 

ID 


(2) 

Best-fit Model 


(3) 
N H 


(4) 

r/fcr/r in 


(5) 
X 2 /DOF 


(6) 

Comments 


7 

13 
19 
21 

26 


P 

BB 
P 
T+P 
BB 


6.0_ .5 
15 3 +1 - 9 

^•"-1.4 

6 4 +07 

—1.0 
AU -0.7 


i AQ +0.09 
1 -° S -0.08 

, -t+0.25 
- 2 - 1 '-0.12 

1 S4+ 03 ' 

*7- = 0.3(t^ 7 /r = 2.52^^ 
, 4^+0:21 

L46 -0.17 


166.0/193 
113.0/108 
17.6/16 
94.4/97 
45.1/52 


BB model also acceptable, but shows systematic residuals 
P model shows somewhat worse residuals 

One-component models are not acceptable, Z = 6.39(> 1.45)Zq 
P and T models also acceptable, but show some systematic residuals 


ID 


Best-fit Model 


A'll 


T/kT/Tin 


cstat/DOF 


Comments 



6 
9 

11* 
20* 
30 
31 

33 



6.75 
5.2; 
6.2 

2.8: 
6.9; 



+0.8 

-1.0 
+0.7 
-0.6 
+3.6 
-1.9 
+4.3 



4.7 



+0.8 
-0.7 



7+0.16 
'-0.16 
-,+D. Id 
z -0.14 
7+0.03 
'-0.01 
-+0.04 
--0.0K 

3 14+1 54 
g+041 
-(>.?(. 

'-0.18 



1.775 

1.825 
0.175 
0.255 



3.795 
2.175 



1375.5/3233f 
1255.7/3233f 
870.6/3233f 
706.2/3233f 
842.8/3233f 
969.8/3233f 
1046.1/3233f 



T model also acceptable, but shows worse residuals 
Soft X-ray spectrum, consistent with star ID 

Soft X-ray spectrum, consistent with faint star ID, Z = 0.285j ] 3 Zq 
T and BB models also acceptable. 
T and BB models also acceptable. 



NOTE. — Col. 1: ID number. Sources denoted by a "*" are not thought to be intrinsic to IC 342. Column 2: Best-fit spectral model from XSPEC for bright sources (> 150 counts). 
For sources with fewer than 400 total counts, the Cash statistic ' Cash 1979) was used instead of \ 2 - All models were fit with a wabs absorption component. A "P" indicates a power-law 
compon ent, a "T" indicates a mekal thermal-plasma component, and a "BB" indicates a multi-color disk component. Note that the Galactic column toward IC 342 is = 3 X 10 21 cm" 2 
IStark et all 19921) . Column 3: Neutral hydrogen absorption column density in units of 10 21 cm" 2 as determined from the best-fit models to the XMM-Newton spectra. Also listed are 
the 90% confidence errors calculated for one parameter of interest (Ax" = 2.7). Column 4: Power-law photon index T, thermal-plasma temperature kT (keV), or multi-color disk inner 
temperature T m as determined from the best-fit models to the XMM-Newton spectra. Also listed are the 90% confidence errors calculated for one parameter of interest (Ax 2 = 2.7). 
Column 5: x" an d degrees of freedom. Sources denoted by a "f" were fit with the Cash statistic on the unbinned data, in which case the "cstaf" value and degrees of freedom are listed. 
Column 6: Comments. 



If X-l truly is in a VH state, and this state has remained 
roughly constant for w 2 yrs, then it marks the longest pe- 
riod over which a VH state has been observed. For compar- 
ison, this would be « 10 x longer than has been observed in 
GX 339-4, a Galactic binary which went into a VH state for 
«3 months (M. Nowak 2003, private communication). Given 
the limited amount of information regarding the VH state in 
general, the longevity of the VH state observed in X-l is plau- 
sibly consistent with current understanding of this state. An- 
other possibility is that both the 2000 ASCA and 2002 XMM- 
Newton observations caught X-l in a VH state, but it varied 
out of the VH state in between. GRS 1915+105, for instance, 
is known to vary in and out of the VH state freq uently a nd has 
a non -negligible duty cycle in the VH state (e.g., Bello ni et alJ 
2000). 

3.3.2. XMMU J034616.0+681115 (IC 342 X-2) 

IC 342 X-2 is the second brightest source in IC 342 and 
also has been extensively studied by ASCA. In contrast to X- 
1, this source was in a low state in 1993 with Fo.5-10 keV ~ 
4.1 x 10" 12 erg cm" 2 s" 1 . Its X-ray spectrum was acceptably 
fit by either an absorbed power-law or diskbb model (098, 
K01). The best-fit absorbed power-law parameters were Nh = 
(14.3 ± 1.6) x 10 21 cm" 2 and a photon index of T = 1.39 ± 
0.10 keV. In 2000, K01 found that the flux of this source had 
nearly doubled to Fo.5-10 keV ~ 7.2 x 10~ 12 erg cm" 2 s" 1 , and 
its spectrum was best fit by an absorbed diskbb model with 
7V H = (18.0 ± 8.0) x 10 21 cm" 2 and T m = 1 .62 ± 0.04 keV. K01 
attributed this to be a spectral transition between a low/hard 
state and a high/soft state. As with X-l, we note that the 
ASCA «3f0 radius beam used to study X-2 would have also 
included source 12 and possibly source 20 (compare with Fig- 
ure lb of K01). The contamination from these sources to the 
ASCA source, however, is only w 3-10% based on the relative 
XMM-Newton source fluxes (and should be even less after the 
strong variability seen from X-2 is accounted for; see below). 

Based on the ASCA spectral fits, we fitted both absorbed 



diskbb and power-law models to the X-ray spectrum of X- 
2. We note that X-2 lies along a chip gap in the MOS2 
detector, resulting in a 50% loss of counts with that de- 
vice. The best-fit parameters for the absorbed diskbb model 
were T m = 2.17: 1 ; 2 ^ keV and 7V H = (15.3+J;|) x 10 21 cm" 2 
(X 2 = 113.0 for 108 degrees of freedom), while those for 
the absorbed power-law model were T = 1.81 ±0.18 and 
7V H = (22.5:^) x 10 21 cm" 2 ( X 2 = 127.3 for 108 degrees of 
freedom). Both models provide acceptable fits to the data, 
although the residuals of the diskbb model exhibit fewer sys- 
tematic deviations (see Figure^. The best-fit parameters are 
different from those obtained during either ASCA observation, 
demonstrating that X-2 has undergone further spectral vari- 
ability. In addition, the XMM-Newton flux has dropped by 
factors of ~1.7 and ^3 compared with the ASCA results of 
1993 and 2000. Thus although this state is spectrally interme- 
diate between the 1993 and 2000 states, it is the lowest flux 
state exhibited by X-2 to date. It is interesting that this new 
low/hard state is better fit by the diskbb model than a power- 
law model, as this indicates that X-2 still has some spectral 
curvature. This curvature is in the opposite sense to that ex- 
pected from a soft thermal component with a hard tail (i.e., 
the standard for X-ray binaries in their high state). The col- 
umn density observed toward this source, however, is a factor 
of «2-3 higher than seen from any of the other sources, sug- 
gesting that X-2 has significant local absorption. It is there- 
fore possible that the apparent curvature is due to complex 
absorption near the source, which is also often detected in lu- 
minous Galactic X-ray binaries. 

3.3.3. XMMU J034648.8+680546 (IC 342 X-3 ) 

IC 342 X-3 is coincident with the nucleus and was clearly 
detected by both the ROSAT HRI and ASCA. The spectral 
properties of X-3 from the 1993 ASCA observation are dis- 
cussed briefly in 098, who found that the source was best 
fitted by an absorbed power law with T = 2.1 ±0.2 and fixed 
7V H = 3.0 x 10 21 cm" 2 . Again, we note that the ASCA tv3'.0 ra- 
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FIG. 4. — X-ra y spectra of the four brightest X-ray sources in IC 342. The upper panel of each plot shows the PN, MOS1, and MOS2 spectra and best-fit 
models (see i|3.3l for details). The lower two panels show the residuals of the fits measured in units of a for the best-fit model (middle) and the next most likely 
model (bottom). 



dius beam used to study X-3 would have also included sources 
19, 24, 26, and 29 (compare with Figure 1 of 098). The con- 
tamination from these sources to the ASCA source is » 44% 
based on the relative XMM-Newton source fluxes. 

We initially fitted both single absorbed power-law and 
thermal-plasma models to the X-ray spectrum of the nucleus. 
Both fits were unacceptable, with the absorbed power-law 
model giving T = 2.67 and N H = 4.0 x 10 21 cm" 2 ( X 2 = 180.4 
for 101 degrees of freedom), and the absorbed mekal model 
giving kT = 2.14 keV, Z = 0.02Z Q , and JV H = 2.4 x 10 21 cm" 2 
(X 2 = 214.9 for 100 degrees of freedom). Neither model ac- 
ceptably fits the continuum shape or the apparent Fe L-shell 
emission complexes visible at «0.8 keV and ssl.O keV (see 
Figure |4}. The presence of Fe L-shell emission at these 
two energies suggests that a large fracti on o f the emission 
is likely to be thermal in origin. In A3. II we were able 
to separate spatially X-3 into point-like and extended com- 
ponents. Thus we should expect that the spectrum of X- 
3 might be complex, consisting of perhaps two or more 



components. We therefore modeled the spectrum with an 
absorbed mefca/+power-law model. This model provides a 
significant improvement over the single-component models 
( x 2 = 94.4 for 96 de grees of freedom), with best-fit values 
of kT = 0.30^$ keV, T = 2.52^;}f, Z = 6.39(> 1.45)Z Q , 
and Nn = (6.8 + [g) x 10 21 cm -2 . An absorbed mekal+mekal 
model with a single metallicity also provides an acceptable 
fit to the data ( X 2 = 97.9 for 96 de grees of freedom), with 
best-fit values of kT x = 0.50^ keV, kT 2 = 4.02+} keV, 
Z = 0.12!°$Z Q , and N H = (5.0! 2 ; 2 ) x 10 21 cm" 2 . However, 
the best-fit abundance for this second model is low com- 
pared to the optically derived value of w 3Z iDiaz & Tosil 
1986; Vila-Costas & Edmunds 1992). Fixing the absorbed 
mekal+mekal model abundance at 3Z Q (or even Z©) results 
in an unacceptable fit with X 2 = 169.9 (x 2 = 120.5) for 97 de- 
grees of freedom. The absorbed mefcaZ+power-law model on 
the other hand provides a consistent measure of the nuclear 
abundance. 

The relative contribution of the mekal component to the to- 
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tal absorbed flux of X-3 is 45% and 1% in the soft and hard 
bands, respectively. The soft-band fraction is in basic agree- 
ment with the physical picture found from our spatial anal- 
ysis ( H3.ll . suggesting that the soft thermal plasma compo- 
nent is likely to be spatially extended and of a starburst ori- 
gin. The complete domination of the power-law component 
above 2 keV, coupled with the fact that we still find signif- 
icant amount extended emission in this band, suggests that 
the hard component is probably complex and may be due to 
an unknown number of circumnuclear X-ray binaries. 6 The 
absorption-corrected 0. 5-10. OkeV luminosity of the extended 
nuclear component is Lx ~ 2.1 x 10 39 erg s -1 ), which is com- 
parable to or larger than our constraint on extended emission 
from the rest of IC 342. 

3.3.4. XMMU J034657. 7+680616 (IC 342 X-4) 

IC 342 X-4 is a bright source located «1.'0 from the nu- 
cleus, which was only resolved by the ROSAT HRI observa- 
tion (ROSAT source 9; BCT93). The source was presumably 
detected by ASCA as well but was blended together with the 
nucleus in those observations, contributing w 31% based on 
the relative XMM-Newton source fluxes. We initially fitted 
X-4 with an absorbed power-law model, obtaining an accept- 
able fit with best-fit parameters of Nh = (4.4+Jq) x 10 21 cm" 2 
and T = (x 2 = 51.8 for 52 degrees of freedom). We 

found that diskbb and mekal models provided equally accept- 
able fits. However, while the diskbb model fit the continuum 
shape best (see the X-4 residuals in Figure 0J, the best-fit ab- 
sorption column value lies just below Galactic. In the case of 
the mekal model, the high temperature would imply that X-4 
must be a young supernova, although it clearly lacks the long- 
term variability or high equivalent- width X-ray emission lines 
characteristic of such sources. Thus it seems that the power- 
law model is the most physically appropriate model. If X-4 is 
a single X-ray binary, then its low photon index may indicate 
that it is in a low/hard state; however, its spectral ambiguity 
leaves room for several other possibilities as well. 

4. CONCLUSIONS 

Our XMM-Newton observation has tripled the number of 
known X-ray point sources in IC 342 and has provided a num- 
ber of new constraints on the X-ray properties of both the 
point-source population and any potential diffuse emission. 

• We find a relatively flat XLF slope of 0.5 ±0.1 for 
IC 342, which is comparable to those found for ac- 
tively star-forming galaxies such as M82 and the An- 
tennae and is inconsistent wit h those of q uiescent spira l 
galaxies fe.g..lKil ffard et al.ll2002HColbert etai]f2 003). 
This implies that IC 342 may have played host to 
a significant amount (> 0.5 M Q yr" 1 ) of recent star- 
formation activity. However, our upper limit on dif- 
fuse X-ray emission associated with hot gas (Lx < 
5.6 x 10 39 erg s" 1 for a 0.5 keV thermal plasma), and 
its infrared luminos ity of L^ M ,„ = 9 x 10 42 erg s -1 
(Shan lev et al.ll200"l ." are more typic al of other Sab 
galax ies such as M58 or M94 (e.g., Eracleous et al. 
2002) rather than luminous starbursts (e.g., Re ad et al. 
1997). Thus IC 342 appears to have the X-ray binary 
population of a starburst, but lacks the hot gas and lumi- 
nous infrared emission that typically accompanies vig- 



orous star formation. This suggests that the mode of 
star formation in IC 342 is probably quite different from 
that of archetypal starburst sources. 

• We find no evidence for short-term variability among 
the X-ray sources in this short observation. We do find 
long-term variability between this observation and the 
1991 ROSAT and 1993/2000 ASCA observations for five 
sources. 

• We find X- 1 to be in a state identical to the one observed 
by ASCA in 2000, after correcting the ASCA measure- 
ment for source contamination. If X-l is in a VH state 
as contended by K02, and it has remained in this state 
between 2000 and 2002, then it would be the longest 
duration VH-state binary ever observed. Alternatively, 
it might be possible that both the 2000 ASCA and 2002 
XMM-Newton observations caught X-l in a VH state 
by chance, but that it has varied out of the VH state in 
between. 

• We find X-2 has decreased in flux by factors of ~1.7 
and ~3 compared to its states observed by ASCA in 
1993 and 2000, respectively, while its spectrum is in- 
termediate between the low/hard 1993 state and the 
high/soft 2000 state. This is the lowest luminosity state 
observed for X-2 to date. 

• We have confirmed that a large fraction («55% and 
«35% in the 0.3-2.0 keV and 2.0-10.0 keV bands, re- 
spectively) of the X-ray emission coincident with nu- 
cleus of IC 342 is associated with an extended, rela- 
tively uniform disk approximately 8" in radius. The 
spectrum of the nucleus is best fit by an absorbed two- 
component model consisting of a thermal plasma with a 
temperature of kT « 0.3 and a power law with a photon 
index of T w 2.53. The relative fluxes of the two spec- 
tral components suggest that the nucleus is complex, 
with a soft extended component contributing «50% of 
the total luminosity. 

Because of their relative distance and luminosity, the X- 
ray binaries in IC 342 provide an excellent laboratory to ad- 
vance our knowledge of ultraluminous X-ray sources even 
further. For instance, a higher signal-to-noise spectrum of X-l 
with XMM-Newton or Constellation-X might provide stronger 
constraints on its nature, especially regarding the origin of 
the spectral softening above 4 keV and an apparent spectrum 
bump around 1.5 keV. Uniform temporal monitoring of X- 
1, X-2, and X-3 could provide more information about the 
frequency and amplitude of spectral transitions, which would 
help place stronger limits on their black hole masses. Finally, 
higher resolution X-ray imaging of the nucleus with Chan- 
dra would allow a decoupling of the extended and point-like 
emission with the two spectral components. 
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sion is from a low-luminosity active nucleus, although there is no evidence 
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